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I . INTRODUCTION 


There  are  a hierarchy  of  models  that  describe  the  phenomenology  of 
double-base  propellants.  Interior  ballistics  codes  accept  propellant 
data  as  inputs,  and  predict  pressure  vs.  time  curves,  muzzle  velocities, 
etc.  Consideration  of  these  models  is  beyond  the  scope  of  this  research, 
and  the  mention  is  made  only  to  set  limits  to  the  capabilities  of  other 
models . 

Among  other  inputs,  propellant  combustion  modelers  try  to  predict 
the  regression  rate  of  a burning  propellant.  A model  to  predict  the 
regression  rate  of  a solid  requires  knowledge  of  many  properties,  such 
as, 


(1)  Solid-state  reactivity  and  heat  release, 

(2)  Gas-phase  chemistry, 

(3)  Gas-phase  transport  properties, 

(4)  Coupling  mechanism  between  the  gas  and 

condensed  phases. 

Little  information  is  available  for  any  of  these  properties,  although 
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recent  work  in  this  laboratory  may  shed  light  on  items  two  and  three  ’ . 
In  particular,  Kooker's  modeling  provides  the  first  realistic  coupling 
of  chemistry  and  hydrodynamics  for  a confined  reacting  gas  mixture. 

The  research  described  in  this  report  was  motivated  by  the  need  for 
pressure  and  temperature  dependent  rate-constant  expressions  for  the 
unimolecular  decomposition  reactions  of  two  free-radical  species,  methoxy 
radical  (CHjO)  and  formyl  radical  (HCQ) . Such  data  are  required  for 

qualitative  kinetic  modeling  of  the  primary  flames  of  simple,  mono- 
nitrate esters.  Models  of  the  primary  flames  for  this  class  of  compounds 
should  provide  transferable  data  for  models  of  the  gas-phase  reactions 
for  burning  propellants. 

Pressure  and  temperature  dependent  rate  constants  for  unimolecular 
reactions  can  be  predicted  by  application  of  the  Rice-Ramsperger-Kassel- 
Marcus  (RRKM)  theory.  The  next  section  reviews  the  major  features  of 
this  theory,  and  subsequent  sections  deal  with  the  application  of  the 
method  to  methoxy  radical  and  formyl  r»diral. 


1.  R.A.  Fifer,  "High  Temperature  P.  •oiysis  of  Methyl  and  Ethyl 
Nitrate",  Seventeenth  Combustion  ■ ■ - titute  Symposium  preprint. 

2.  D.A.  Kooker,  "Transient  Laminar  Fi.  :>o  Propagation  in  u-n ,<*d 
Premixed  Gases;  Numerical  Prediction.",  Seventeenth  Combustion 
Institute  Symposium,  preprint. 


5 


II.  REVIEW  OF  THEORY 


A recent  report  reviews  the  development  and  application  of  RRKNF*. 
More  complete  discussions  of  the  theory  for  unimolecular  reactions  can 
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be  found  in  texts  by  Forst  and  Robinson  and  Holbrook  . Several  points 
about  the  theory  require  that  a modest  overview  precede  discussing  the 
results  of  this  project. 

The  Lindeman  theory  forms  the  basis  of  almost  all  modern  treatment 
of  unimolecular  reactions,  The  mechanism  consists  of  three  steps: 

(1)  An  energizing  collision  between  the  reactant,  A,  and 

a collision  partner,  M,  produces  an  energized  molecule 
* 

‘ A . 

* 

A+M  — ar«  A +M.  (1) 

The  energized  molecule  possesses  energy  greater 
than  a critical  energy,  Eo*  The  critical  energy 

is  the  minimum  energy  the  species  A must  possess 
in  order  to  undergo  a chemical  transformation. 

(2)  Energized  molecules  are  de-energized  by  collision. 

* ^2 

A +M  ^ A+M  (2) 

The  Lindeman  theory  assumes  that  every  collision 
* 

of  A leads  to  de-energizing. 


(3)  A time-lag  exists  between  the  energizing  of  A and 
the  unimolecular  reaction  of  the  energized  molecule. 

* ^3 

A _ ►B+C.  (3) 

The  chemical  rate  coefficients  are  assumed  to  be  constants.  By  applies 

★ 

tion  of  the  steady-state  hypothesis  to  the  concentration  of  A , an 
equation  for  the  overall  rate  of  reaction  can  be  written, 

v - k3(A*)  = k1k3(A)  (MVCk2(M)  -k3)  . (4) 

T.  G.F.  Adams,  ''An  Analysis  of  the  Pressure  Dependence  of  Nitrate 
Ester  Thermal  Decomposition",  ARBRL-TR-02106.  (AD  #A063513) 

4.  W.  Forst,  "Theory  of  Unimolecular  Reactions",  Academic  Press, 

New  York,  1973. 

P.J.  Robinson  and  K.A.  Holbrook,  "Unimolecular  Reactions",  Wiley 
Interscience,  New  York,  1972. 
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There  are  two  important  limiting  cases.  At  high  pressures,  k,(M)  » k_, 
so  that 


V.  3 (A)  3 k«  <A)  ‘ (5) 

At  low  pressures,  k3  » k2  (M) , and 

Vbim  3 ki (A) (M) , (6) 

a rate  that  depends  on  the  second-order  energy-transfer  rate  constant. 

Lindeman's  theory  predicts  that  the  unimolecular  rate  coefficient 
is  not  constant  for  low  concentrations,  but  is  constant  at  high-pres3ures. 
These  predictions  are  qualitatively  correct,  but  quantitative  predictions 
based  on  this  theory  are  never  correct,  principally  because  the  rate 
coefficients  kj,  k2,  k3  are  not  constants.  Modifications  to  Lindeman's 

theory  led  to  development  of  the  RRKM  theory.  The  predictions  of  RRKM 
theory  agree  with  experiment.  A review  of  the  application  of  the  RRKM 
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method  has  been  given  . The  emphasis  here  is  that  the  reactions  proceed 
with  variable  molecularity  - bimolecular  at  low  pressure,  unimolecular 
at  high  pressure,  and,  in  the  transition  region  between  these  two  limits, 
with  a variable  molecularity.  Figure  1 illustrates  the  various  regions, 
showing  a plot  of  the  unimolecular  rate  coefficient,  k(m^,  versus  P for 

a hypothetical  reaction.  It  is  the  purpose  of  this  report  to  describe 
the  pressure-dependence  of  the  decomposition  reactions  for  the  methoxy 
and  formyl  radicals,  and  to  provide,  in  the  case  of  methoxy,  a pressure- 
dependent  rate  expression  for  the  transition  region. 

The  following  sections  describe  the  methoxy  radical  decomposition 
and  the  formyl  radical  decomposition.  To  predict  the  rate  constants 
for  these  reactions,  it  is  necessary  to  have  information  concerning 
the  molecular  properties  of  the  reacting  system.  Specifically,  one  needs 
to  know: 

1.  The  critical  energy,  E , for  the  decompose  Lon 
reaction, 

2.  The  vibrational  frequencies  for  the  reactant 
molecule  and  the  energized  molecule, 

3.  The  molecular  structure  parameters,  bond  lengths 
and  angles,  for  the  reactant  molecule  and  the 
energized  molecule. 


III.  A PREDICTED  RATE  CONSTANT  FOR 
METHOXY  RADICAL  DECOMPOSITION 

The  methoxy  radical,  CHjQ,  is  a rara  avis.  Attempted  spectroscopic 

studies  of  the  molecule  have  produced  a spectral  signature,  but  not 

enough  data  to  provide  a molecular  structure  for  the  compound  A 
theoretical  quantum  chemical  study  of  the  electronic  structure  of  the 

radical  has  provided  an  estimate  of  the  bond  lengths  and  angles7. 
Yarkony,  et  al.,  found  that  the  bond  lengths  for  CHjO  were  approximately 

equal  to  the  respective  bond  lengths  in  methanol.  Further,  the  theo- 
retical results  imply  that  the  C-0  bond  is  not  perpend.'  .alar  to  the 
plane  of  the  three  hydrogens,  but  is  b'ent  three  degrees  toward  the  two 
equivalent  hydrogens.  This  staggered  arrangement  breaks  the  C3  symmetry 

of  the  methyl  group,  leaving  only  a plane  of  symmetry.  The  predicted 
structure  provides  the  data  needed  to  compute  the  moments  of  inertia 
for  the  reactant.  In  addition,  since  the  predicted  structure  is  so 
close  to  that  of  methanol,  .the  vibrational  frequencies  for  the  reactant 
are  assumed  equal  to  those  for  methanol.  These  estimates  and  the  com- 
puted moments  are  presented  in  Table  1. 

Because  the  theoretical  study  of  methoxy  included  no  details  of 
the  potential  energy  hypersurface,  no  structure  was  predicted  for  the 
transition  state.  Decomposing  methoxy  to  form  formaldehyde  and  a 
hydrogen  atom  requires  several  structural  changes: 

1.  The  stretching  of  a C-H  bond  to  an  effectively 
infinite  length, 

2.  The  flattening  of  the  methyl  group,  leaving  a 
planar  CH2  terminus, 

3.  The  shortening  of  the  C-0  bond. 

None  of  these  changes  will  be  complete  at  the  transition  state,  however. 
The  model  used  here  is  obtained  by  making  the  following  assumptions. 


6.  H.E.  Radford  and  D.K.  Russell,  "Spectroscopic  Detection  of  Methoxy 
(CH-O)",  J.  Chem.  Phys.  66,  2222  (1977). 

7.  D.R7  Yarkony,  H.F.  Schae??er,  III  and  S.  Rothenberg,  "Geometries 
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of  the  Methoxy  Radical  (X  E and  A A,  States)  and  the  Methoxide  Ion", 
J.  Am.  Chem.  Soc. , 96,  656  (1974). 
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1.  The  ”in-plane"  C-H  bond  is  stretched  to  a bond 
length  of  3.S  A. 

2.  The  CO  bond  length  is  shortened  to  1.3  A,  between 
the  values  for  methoxy  and  formaldehyde. 

3.  The  O-CH^  bond  angle  is  reduced  by  approximately 
ten  degrees. 

Table  1.  Molecular  Data  for  Methoxy  Radical  and  Excited  Complex8. 


Property 

Radical 

Complex 

RCH 

1.08 

3.5 

RC0 

1.44 

1.3 

0OCH 

109° 

109° 

E0 

24  kcal  mol  * 

rxx 

19.5 

25.3 

! I 

19.53 

17.1 

1 yy 

i *zz 

3.62 

11.3 

1 

2977 

2980 

*2 

2977 

2840 

V3 

2844 

0000 

V4 

1477 

1400 

VS 

1455 

1000 

V6 

1430 

300 

V7 

1171 

300 

V8 

1112 

850 

V9 

1029 

850 

g 

Bond  distances  are  given 
vibrational  frequencies 

in  A,  moments  of  inertia 
in  wave  numbers. 

in  amu- 

Stretching  of  the  CH  bond  causes  one  of  the  C-H  stretching  frequencies 
to  approach  zero,  and  reduces  four  of  the  six  methyl  bending  modes. 
Estimates  of  the  vibrational  frequencies  for  the  energized  molecule  and 
the  computed  moments  of  inertia  are  given  in  Table  1. 


The  vibrational  frequencies  and  moments-of-inertia  assumed  for  the 
methoxy  radical  and  the  transition  species  are  used  to  compute  the 
rotational  and  vibrational  partition  functions.  In  addition,  a value 
for  the  high-pressure  Arrhenius  factor  can  be  estimated. 

To  perform  the  RRKM  integration,  these  factors  must  be  supplemented 
by  a value  for  the  energy  barrier  to  the  unimolecular  reaction.  Figure 
2 shows  a schematic  plot  of  potential  energy  versus  an  arbitrary  reaction 
coordinate.  The  potential  energy  barrier,  indicated  by  E,  is  the  simple 
difference  in  the  electronic  energy  for  the  two  configurations.  This 
number  must  be  corrected  to  include  the  difference  in  the  zero  point 
vibrational  energies  for  the  two  structures,  giving  a value  for  the 
critical  energy,  Eq.  Integration  of  the  RRKM  equation  for  the  unimolec- 
ular rate  constant  requires  a knowledge  of  the  critical  energy.  Since 
no  experimental  value  for  this  parameter  is  available,  an  estimate  must 
be  made. 

Among  thermochemi cal  data  available  for  methoxy  is  the  dis- 

8 -1 
sociation  energy  relative  to  formaldehyde  , D(H-CH20)  * 29  kcal  mol'  . 

It  is  common  to  find  empirical  estimates  of  activation  energies  that  are 

lower  than  the  dissociation  energy,  often  by  five  kcal  mol  , or  more. 

Applying  this  empirical  rule,  a value  of  24  kcal  mol-1  is  assigned  to  the 
activation  energy,  E , which  is  assumed  to  equal  the  critical  energy, 

a 

A thermochemical  argument  by  Batt  and  McCulloch  estimates  E = 25  kcal 

- 1 a 

mol  . The  uncertainty  associated  with  the  two  methods  of  estimation 

suggest  that  this  near  agreement  is  fortuitous.  Predictions  of  kyni 
for  the  methoxy  radical  decomposition  are  made  by  assuming  Eq  = 24  kcal 
mol"1.  An  alternative  theoretical  method  for  estimating  Eq  is  pre- 
sented in  the  discussion  of  the  results  for  the  formyl  radical  decom- 
position. 

To  complete  describing  the  analytical  method  used  to  predict  the 
-'ate  constant,  it  is  noted  that  the. sums  and  densities  of  quantum  g 

states  were  calculated  via  the  Whitten-Rabinovitch  semiclassical  method  . 
A discussion  of  the  application  of  the  method  is  given  in  references 
(3)  and  (5). 


J.  Heicklen,  "Gas  Phase  Reactions  of  Alkylperoxy  and  Alkoxy 
Radicals",  Adv.  Chem.  Ser.,  No.  76,  23  (1974). 

,’.Z.  Whitten  and  B.S.  Rabinovitch,  "Accurate  and  Facile  Approxi- 
mation for  Vibrational  Energy-Level  Sums",  J.  Chem.  Phys.,  38,  2466 

( '!  C|fi  ; ) . 
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Figure  2.  A schematic  representation  of  the  variation  of  electronic 
energy  with  reaction  progress.  The  critical  energy  is 
shown  as  the  difference  between  the  lowest  total  energy 
state  of  the  reactant  and  the  transition  state. 


The  rate-constant  for  the  methoxy  radical  decomposition  was  com- 
puted at  four  temperatures:  900*,  1000°,  1100®,  and  1200°K.  The 
results  of  the  four  calculations  are  shown  ir.  l.ole  2.  The  kinetics 
data  is  summarized  in  figure  3,  where  the  hatched  area  identifies  the 
pressure  range  of  interest  for  the  pyrolysis  model.  This  area  obviously 
lies  in  the  fall -off  portion  of  the  rate-constant  curves;  thus  neither 
the  high-pressure  nor  the  low-pressure  limit  for  the  rate  expression 
will  he  valid  for  the  modeling. 

The  pressure  dependence  of  the  unimolecular  rate  constant  demon- 
strates the  competition  between  the  collisional  energy  transfer  and 
intramolecular  relaxation.  These  two  processes,  essentially  decoupled 
i.i  the  low-  and  high-pressure  limits,  overlap  in  the  intermediate  part 
of  the  "fall-off"  region.  Troe  has  suggested  that  the  difficult  analy- 
sis of  the  intermediate  fall-off  curve  shoul^he  based  on  information 
derived  from  the  two  limiting  rate  constants1  . Troe  derived  a general 
first-order  rate  constant  k/k  as  a measure  of  pressure.  His  analysis 
concluded  that  the  reduced  fall-off  curves  do  not  depend  strongly  on 
the  theoretical  model  of  the  reaction.  This  suggests  that  even  simple 
theoretical  models  are  sufficient  to  predict  reduced  fall-off  curves 
useful  for  extrapolating  experimental  data, 

Troe's  objective  was  to  predict  the  entire  fall-off  region  well 
enough  to  extrapolate  data  from  one  limiting  case  to  another.  The 
goal  of  the  present  research  is  to  derive  a rate  constant  expression 
that  adequately  describes  a relatively  small  part  of  the  rate  constant 
curve.  1'ollowing  Troe,  a reduced  rate  constant  k/km  is  defined. 

Assuming  that  can  be  expressed  by  an  Arrhenius-expression 


k = A e 

CO  CO 


VI,T 


(7) 


there  are  two  alternatives  for  defining  A^: 

1.  Let  Am  = (T)  for  the  temperature 

of  interest. 

2.  Define  a single,  average  high-pressure 
frequency  factor  = < Aw(T)  >.j,. 


10.  J.  Troe,  "Fall-Off  Curves  of  Unimolecular  Reactions",  Der.  Bunsenges. 
physik.  Chem.  , 78,  474  (1974). 
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Table  2.  Predicted  k . for  Methoxy  Uadicala. 

uni 


Pressure 

(torr) /Temperature 

900 

(°K) 

1000 

1100 

1200 

10 

4.5(5) 

1.5(6) 

3.9(6) 

8.3(6) 

100 

4.4(6) 

1.5(7) 

3.8(7) 

8.2(7) 

760 

2.9(7) 

9.9(7) 

2.7(8) 

5.9(8) 

1000 

3.6(7) 

1.3(8) 

3.4(8) 

7.6(8) 

1520 

5.1(7) 

1.8(8) 

4.9(8) 

1.1(9) 

3800 

9.6(7) 

3.5(8) 

1.0(9) 

2.4(9) 

5000 

1.1(8) 

4.2(8) 

1.2(9) 

2.9(9) 

7600 

1.4(8) 

5.4(8) 

1-6(9) 

3.9(9) 

10000 

1.6(8) 

6.3(8) 

1.9(9) 

4.7(9) 

11400 

1.7(8) 

6.7(8) 

2.1(9) 

5.1(9) 

15200 

1.9(8) 

7.7(8) 

2.4(9) 

6.0(9) 

20000 

2.1(8) 

8.6(8) 

7.8(9) 

7.0(9) 

30000 

2.4(8) 

1.0(9) 

3.3(9) 

8.5(9) 

38000 

2.6(8) 

1.1(9) 

3.6(9) 

9.4(9) 

50000 

2.7(8) 

1.2(9) 

4.0(9) 

1.0(10) 

76000 

3.0(8) 

1.3(9) 

4.5(9) 

1.2(10) 

100000 

3.1(8) 

1.4(9) 

4.9(9) 

1.3(10) 

114000 

3.2(8) 

1.4(9) 

5.0(9) 

1.4(10) 

150000 

3.3(8) 

1.5(9) 

5.3(9) 

1.5(10) 

200000 

3.4(8) 

1.5(9) 

5.6(9) 

1.6(10) 

*1 

14875 

16122 

18600 

21194 

28012 

32808 

37850 

43128 

^2 

2.466 

3.09 

3.9 

4.95 

log 

15.39 

15.5 

15.56 

15.61 

a Notation:  4.1(5)  implies 

4.1  X 105. 

Statistical 

mechanical  parameters 

are  dimensionless. 


Figure  3.  Variation  of  kun.  with  pressure  for  the  decomposition  of 


methoxy  at  three  temperatures,  1000DK,  1100*K  and  1200®K. 
The  hatched  area  identifies  the  physical  conditions 
encountered  in  Fifer's  experiments  (1). 
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The  energy  in  the  Arrhenius  expression,  Equation  7,  is  assumed  to  be  the 
critical  energy,  Eq.  Table  3 contains  the  results  of  computing  k/k^  for 

each  choice  of  A^,  at  various  pressures  and  temperatures  corresponding  to 

those  encountered  in  the  pyrolysis  experiments.  Figure  4 presents  a 
plot  of  log(k/kw)  versus  log  P for  the  case  A^  = < Aoo(T)>.  The  plot 

shows  that,  to  reasonable  accuracy,  a single  line  could  be  used  to 
describe  these  data.  Table  4 lists  the  average  k/k^  ratios  for  pressures 

in  the  range,  2 atm  < P < 20  atm.  A plot  of  log  (k/km)  versus  log  (P) , 

using  the  average  values  for  the  reduced  rate  constant  has  a slope,  m, 
of  0.65.  To  derive  the  final  expression  for  the  unimolecular  rate 
constant,  the  following  equation  must  be  integrated. 

d log  ■ m d log  P.  (8) 


The  integration  has  a lower- limit  at  P^  ■ 1520  torr,  where  (k/k^)^  ■ 


.0095.  The  integrated  expression  is 

logfk/k^)  - logfk/k^  - m(log  P - log  P^  . (9) 

Rearranging, 

k/kM  a ft/kJi  (P/P±)m  (10) 

- .0095  (P/1520)0,65,  (11) 

so  that 

k - ^(.0095)  (P/1520)0,65,  (12) 

where  km  ■ 3.31  * 1015  exp  (-24000/RT).  (13) 


Kineticists  often  prefer  that  the  rate-constant  expression  replace  the 
pressure  variable  with  a density  variable.  Performing  this  manipula- 
tion yields  the  following  expressions  for  the  rate  constant, 

0.65 

702-44  ) °4) 

In  order  to  evaluate  the  validity  of  the  derived  rate-constant 
expression,  values  of  kuni  predicted  using  this  expression  are  compared 

with  the  results  of  the  original  RRKM  integration.  The  data  presented 
in  Table  5 show  that  the  derived  rate-constant  expression  is  an  excellent 
approximation  for  describing  the  unimolecular  decomposition  rate  for 
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1 


Table  3.  (kuni/kj 

for  Two 

Choices  of  A . 

op 

I.  Pre-exponential  factors 

taken  from  Table 

2. 

Temperature 

900 

1000 

1100 

1200 

Pressure  (torr) 

1520 

.013 

.01 

.0081 

.0064 

3800 

.025 

.02 

.017 

.014 

5000 

.03 

.024 

.02 

.017 

7600 

.037 

.031 

.027 

.023 

10000 

.042 

.036 

.032 

.027 

11400 

.045 

.038 

.034 

.030 

15200 

.05 

.044 

.04 

.035 

I I . Pre-exponential 

factor 

assumed  to  be  the 

average 

of  the  values 

from  Table  2. 

900 

1000 

1100 

1200 

1000 

.0077 

.0068 

.0062 

.0055 

1520 

.011 

.0097 

.009 

.008 

3800 

.02 

.019 

.018 

.017 

5000 

.024 

.023 

.022 

.021 

7600 

.030 

.029 

.029 

.028 

10000 

.034 

.034 

.035 

.034 

11400 

.036 

.036 

.038 

.037 

15200 

.041 

.042 

.044 

.043 

Table  4.  Variation  of  Reduced  Rate  Constant  With  Pressure  for  CH^O. 


ssure  (torr) 

<kuni/k. 

1520 

.0095 

3800 

.0185 

5000 

.0225 

7600 

.029 

10000 

.034 

11400 

.038 

15200 

.0425 
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Table  5.  Comparison  of  RRKM  Predictions  for  k , 

uni 

Reduced  Pull -off  Equation. 


with  Predictions  of 


1100 

Formula 

4.9E8 

S.3E8 

1.0E9 

1.1E9 

1.6E9 

1.5E9 

2.1E9 

2.0E9 

2.4E9 

2.4E9 

the  methoxy  radical  within  the  pressure  and  temperature  ranges  specified 
by  the  pyrolysis  experiment.  In  addition,  the  simple  form  of  the 
derived  expression  recommends  its  use  in  kinetic  models. 

IV.  A PREDICTED  RATE  CONSTANT  FOR  FORMYL  RADICAL  DECOMPOSITION 
Before  discussing  the  prediction  of  kun^  for  the  reaction, 

HCO  -*■  H + CO, 

it  should  be  noted  that  the  formyl  radical  may  not  satisfy  the  criterion 
for  applying  RRKM  theory,  The  breaking  of  a carbon-hydrogen  bond  in 
molecules  with  atoms  other  than  just  carbon  and  hydrogen,  is  often  found 
to  be  slower  than  anticipated.  This  slow-down  in  the  specific  rate  of 
reaction  is  rationalized  as  an  effect  of  inefficient  intramolecular 
transfer  of  vibrational  energy,  a manifestation  of  the  high  frequency 

for  the  C-H  stretching  motion11.  Since  no  dynamic  study  has  yet  demon- 
strated that  the  molecule  is  a non-RRKM  molecule,  the  RRKM  method  will 
be  applied  to  analyze  the  decomposition  reaction. 

lTT  W.C.  Ha  so’,”  private  communication. 


In  contradistinction  to  the  methoxy  radical,  there  is  excellent 
experimental  data  characterizing  the  formyl  radical.  Both  a molecular 
structure  and  a complete  vibrational  spectrum  are  known  for  the  formyl 
radical  ground  electronic  state.  Furthermore,  a theoretical  investi- 
gation of  the  electronic  structure  of  this  radical  includes  a self- 
consistent  field  estimate  of  the  potential-energy  hypersurface  for  the 
decomposition  reaction.  This  surface  provides  a molecular  structure 
for  the  activated  complex,  as  well  as  an  estimate  for  the  critical 
12 

energy  . Given  the  structure  of  the  complex,  it  is  an  easy  matter  to 
estimate  reasonable  values  for  the  vibrational  frequencies  of  the 
complex.  The  data  used  in  this  study  are  presented  in  Table  S. 

A thorough  theoretical  study  of  the  formyl  radical  potential 
energy  hypersurface  predicts  a critical  energy  lower  than  that  predicted 

by  Bruna,  et  al.  The  data  from  the  more  recent  theoretical  work  will  be 

* 

used  in  studies  of  the  formyl  radical  decomposition  reaction  . 

The  RRKM  calculations  were  performed  with  the  same  techniques  used 
for  the  methoxy  decomposition.  Five  different  integrations  were  made, 
corresponding  to  five  choices  of  temperature:  900* , 975® , 1050', 

1125*,  and  1200°K.  Statistical  mechanical  parameters  for  the  calcu- 
lations are  given  in  Table  6,  while  the  predicted  rate  constants  for 
various  temperatures  and  pressures  are  given  in  Table  7.  These  data 
are  displayed  in  Figure  5.  The  reaction  is  controlled  by  bimolecular 
collision  processes  for  all  realistic  values  of  pressure  for  the  temp- 
erature range  considered.  An  expression  for  the  specific  rate  for  the 
bimolecular  process  was  given  in  equation  (6), 

vblm  - kjfAXM)  . (6) 

To  determine  an  expression  for  kj,  the  RRKM  results  for  I1C0  were 

oxtended  to  a pressure  of  1 torr,  giving  an  approximate  variation  of 

with  temperature.  These  data  are  given  in  Table  8.  Using  the  value 
of  k1  for  T-10506,  a pre-exponential  factor  can  bo  derived  which  can 

be  expected  to  give  good  predictions  for  the  othor  values  of  kj.  If  It 

is  assumed  that  the  enerpv-  ln  iiie  exponent  is  just  the  critical  energy, 
the  Arrhenius  u.,1:r'v,t,ioa  :ui>;  the  form 

lTi  F7j7  Bruna,  S.  Peyerimhoff,  and  R.J.  Buenker,  "Ab  initio  Study  of 
the  Structure,  Isomers  and  Vertical  Electronic  Spectrum  of  the 
lormyl  Radical  HCO",  J.  Mol.  Structure,  32,  217  (1976). 

* See  G.F.  Adams,  G.D.  Purvis  and  R.J.  BartTett,  "Double-Excitation 
Many-Body  Perturbation  Theory  Calculations  of  Potential  Energy 
Surfaces.  I.  Formyl  Radical",  AR-BRL  Technical  Report,  in  press. 


Table  6.  Molecular  Data  for  Formyl  Radical  and  Excited  Complex  . 


Property 


Radical 

1.079 


Complex 

1.607 


1.19 

1.13 

130.00 

18.0  kcal  mol"1 

130.00 

11.75 

12.65 

10.86 

10.58 

0.89 

2.06 

2700.00 

1820.00 

2150.00 

1080.00 

100.00 

aBond  distances  are  given  in  angstroms,  momonts  of  inertia  are  given  in 
2 

amu-A  , and  vibrational  frequencies  are  given  in  wave  numbers. 

Table  7.  Statistical  Mechanical  Parameters  for  Formyl  Radical. 


log  A 


900 

975 

1050 

1175 

1200 

4270 

4814 

5380 

5966 

6573 

6646 

7494 

8375 

9288 

10232 

1.3 

1.37 

1.45 

1.53 

1.62 

14.82 

14.84 

14.87 

14.89 

14.90 

k.  - A erp('-E  /RT) 

4 V.'  V 

■ S.SE8  exp f-1 8000/ 1.987T)  . (15) 

Vulues  of  k.  predicted  by  this  e,\p ion  for  et).ch  of  tho  temperatures 
considered  are  given  in  Table  9.  The  agreement,  with  tho  extrapolated 
low-pressure  rate  constants  is  good.  Given  the  . rpress^vj  (15)  foi  k., 
tho  unimolecular  rate  constant  is 

kuni"MP)-  (16> 

This  expression  accurately  reproduces  the  RRKM  predictions  of  tho  uni- 
mplocular  rate  constant. 


••  •!* 


HCO 


T *12  00* 


l;  inure  5.  Variation  of  k ^ with  pressure  for  the  decomposition  of  formyl 

radical  at  three  temperatures,  900#K,  1 050* K , and  1200°K. 

The  hatched  area  identifies  tho  physical  conditions 
encountered  in  Fifor's  experiments  (1), 


Table  8.  Predicted  kuni  for  Formyl  Radical.® 
Pressure  (torr) /temperature  ('K) 


900 

975 

1000 

2.4(7) 

5.2(7) 

sooo 

1.2(8) 

2.5(8) 

10000 

2.3(8) 

4.9(8) 

20000 

4.4(8) 

9.5(8) 

50000 

9.6(8) 

2.1(9) 

75000 

1.3(9) 

2.9(9) 

100000 

1.6(9) 

3.6(9) 

150000 

2.1(9) 

4.7(9) 

250000 

2.7(9) 

6.2(9) 

500000 

3.6(9) 

8.4(9) 

750000 

4.1(9) 

9.6(9) 

106 

4.4(9) 

1.0(10) 

1.5X10* 

4.7(9) 

1.1(10) 

4 Xl0b 

5.2(9) 

1.3(10) 

107 

5.5(9) 

1.3(10) 

1050 

1125 

1200 

9.9(7) 

1.7(8) 

2.8(8) 

4.9(8) 

8.5(8) 

1.4(9) 

9.5(8) 

1.7(9) 

2.7(9) 

1.8(9) 

3.2(9) 

5.3(9) 

4.1(9) 

7.3(9) 

1.2(10) 

5.7(9) 

1.0(10) 

1.7(10) 

7.0(9) 

1.3(10) 

2.1(10) 

9.3(9) 

1.7(10) 

2.8(10) 

1.3(10) 

2.3(10) 

3.9(10) 

1.7(10) 

3.2(10) 

5.6(10) 

2.0(10) 

3.7(10) 

6.5(10) 

2.2(10) 

4.1(10) 

7.2(10) 

2.4(10) 

4.5(10) 

7.9(10) 

2.7(10) 

5.2(10) 

9.3(10) 

2.8(10) 

5.5(10) 

9.9(10) 

dotation:  2.4(7)  implies  2.4X107. 


Table  9.  Low  Pressure  Rate  Constants  for  Formyl  Radical. 


Temperature 

900 

975 

1050 

~rns 

—nm — 

Extrapolated  RRKM 

2.4(4) 

5.15(4) 

9.86(4) 

1.72(5) 

2.8(5) 

Calculated  From 

liquation  (15) 

2.34(4) 

5.08(4) 

9.86(4) 

1.75(5) 

2.9(5) 

Fifor  reanalyzed  the  RRKM  results  replacing  pressure  with  density. 
For  a given  density,  ho  estimates  an  effective  activation  energy  for 
the  bimolecular  reaction,  and  chooses  an  appropriate  pre-exponential 
fuctor.  His  analysis  gives  log  A of  17.0  and  activation  energy  of 

20.1  kcul  mol'1.  He  then  calculates  log  k100Q  ■ 12.6,  which  is  greater 

than  u value  that  can  be  obtained  from  published  data  on  the  revorse 
1 3 

reaction  , log  k1QOO  » 11.1. 

TT  D.L.  Baulch,  D.O.  Orysdale,  J.  Duxbury  and  S.J.  Grant,  "Evaluated 
Kinetic  Data  for  High  Temperature  Reactions,  Vol.  3",  Butterworths, 
London  (1976). 


Most  of  the  difference  between  the  two  estimates  can  be  removed 
by  considering  the  present  treatment  of  collisional  activatio  Two 
variables  that  affect  the  collisional  activation  are  the  collision 
diameter,  a,  and  the  collisional  efficiency,  X.  The  collision  diameter, 
o,  was  assumed  to  be  10  angstroms  for  both  the  methoxy  and  formyl 
decompositions.  This  value  is  probably  too  high  for  the  methoxy  study, 
and  it  is  certainly  too  high  for  the  formyl  radical  investigation. 
Decreasing  <3  to  five  angstroms  would  reduce  the  predicted  rate  constants 
by  a factor  of  4, 

Failure  of  the  strong-collision  assumption  damages  these  efforts  to 
predict  the  rate  constant  for  the  decomposition  of  formyl  radical.  Simply 

stated,  strong  collisions  are  those  which  excite  a molecule  from  a reactant 

* 

state,  A^,  to  a state  of  the  energized  species,  A^,  from  which  reaction  may 

occur.  In  present  terms  the  strong  collision  assumption  is  X»l.  The 
necessity  to  abandon  the  strong  collision  hypothesis  has  been  realized  only 
since  the  1960's,  and  theoreticians  have  concerned  themselves  with  studies 
of  weak  collision  effects.  These  efforts  have  led  to  recent  research  by 

Troe^4*^  and  Stace  and  Murrell^. 

Stace  and  Murrell  developed  an  ingenious  model  to  study  the  effects  of 
pressure  and  bath  gas  on  the  collisional  excitation  efficiency  for  small 
molecules  capable  of  unlmolecular  decomposition.  The  results  of  their 
research  will  be  discussed  in  the  next  report  of  this  series.  Here  the 
consequences  of  Troe's  efforts  are  considered. 

Troe  calculates  the  rate  constants  for  thermal  unimolecular  reactions 
in  the  low  pressure  limit  via  the  solution  of  the  master  equation.,  These 
rate  constants,  k , are  expressed  in  terms  of  a strong  collision  rate 

S C 

constant,  k^  , and  a weak  collision  efficiency,  6^; 


Troc  transforms  this  expression  into  a series  of  complex  equations.  The 
collisional  efficiency  is  the  product  of  a total  cross-sectional  factor, 

6 , such  that, 

E 

IT.  ,J.  Troc,  "Theory  of  Unimolecular  Reactions  at  Low  Pressures.  1. 

Solutions  of  the  Master  Equation.",  J.  Chem.  Phys.,  66,  474S  (1977). 
15.  J.  Troc,  "Theory  of  Unimolecular  Reactions  at  Low  Pressures,  II. 

Strong  Collision  Rate  Constants.  Applications",  J.  Chem.  Phys. 

I(>.  A., I.  Stace  and  J.N.  Murrell,  "Molecular  Dynamics  and  Chemical 

Reactivity.  A Computer  Study  of  Iodine  Atom  Recombination  Under 
High  Pressure  Conditions",  Molecular  Physics,  33,  1,  (1977). 


* ftriTfr  yp«T  - 


(18) 


' Ur  ' 

c at 


The  cross-sectional  factor  is  given  by  a ratio  of  an  energy  dependent 
collision  number  and  a hard-sphere  collision  factor,  while  the  step- 
wise factor  is  taken  to  be 


ec  = - < AE  >/kTF (E) , 


where  < AE  > is. the  average  internal  energy  transferred  per  collision, 
and  F(E)  is  a term  that  corrects  for  the  energy  dependence  of  the 
density  of  states.  In  addition  to  these  relatively  simple  expressions 
for  the  collision-factors,  Troe  developed  a strong-collision  rate 
constant  expression  that  includes  factors  to  correct  for  changes  in 
the  density  of  states  as  a function  of  temperature.  Troe's  important 
result,  however,  is  his  demonstration  that  the  overall  collision 
efficiency,  6 is  a strong  function  of  temperature.  He  finds  that 


at  temperatures  near  3Q0°K,  collisional  efficiencies  have  values  in 
the  range  0.3-0. 5,  a result  which  agrees  with  the  available  experimental 
data.  At  1000°K,  the  collisional  efficiency  has  been  reduced  to  a value 
near  0.1,  and  corresponding  decreases  occur  at  higher  temperatures. 


Troe's  work  demonstrates  that  high  temperature  studies  of  low 
pressure  dissociation  rate  constants  correspond  to  the  weak  collision 
limit  of  energy  transfer,  where  the  "strong  collision  assumption"  of 
unimolecular  rate  theory  has  to  be  abandoned.  The  experimental  con- 
ditions of  formyl  radical  decomposition  that  interest  us  are  high  tem- 
perature, modest  pressure.  Consequently,  it  is  expected  that  the  strong 
collision  assumption  is  not  valid  for  the  experimental  conditions,  nor 
is  the  collisional  efficiency  as  low  as  Troe's  estimates  suggest.  A 


reasonable  guess  would  assign  a value  to  8 of  0.2-0. 3. 

c 


The  combination  of  reduced  collision  diameter  and  reduced  collisional 
efficiency  will  reduce  the  predicted  formyl  rate  constants  by  approxi- 
mately an  order-of-magnitude . It  is  not  fruitful  to  attempt  a better 
cstimute  of  the  effect  since  Troe  studied  only  temperature  effects.  A 
more  complete  treatment  of  the  formyl  radical  decomposition  will  be 
presented  in  the  next  report  in  this  series,  wherein  both  temperature 
and  pressure  effects  will  be  re-evaluated. 


SUMMARY 


This  is  a first  effort  to  predict  the  rate  constants  for  uni- 
molecular decomposition  reactions  of  small  molecular  species.  The  en- 
couraging results  suggest  that  the  techniques  applied  ure  useful  for 
medium- sized  molecules,  such  a3  the  methoxy  radical.  However,  it  seems 


T 1 — ?-*- 
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clear  that  simple,  direct  application  of  the  RRKM  method  to  small, 
three-  and  four-atom  systems  will  give  results  that  are  of  order-of- 
magnitude  accuracy  at  best.  The  next  report  in  this  series  will  present 
two  additional  analyses  of  the  formyl  radical  decomposition,  with  the 
hope  that  more  accurate  results  will  be  obtained, 

The  calculations  reported  assumed  that  the  diluent  gas  was  a strong 
collider.  Before  applying  the  data  in  this  report,  a modeler  should 
ascertain  that  the  strong-collision  assumption  applies  to  the  experi- 
mental process  studied.  For  those  cases  where  the  strong-collision 
assumption  fails,  linear  scaling  of  the  data  contained  here  should 
provide  a reasonable  estimate  for  the  rate  constant. 

The  greatest  uncertainty  in  these  calculations  is  the  critical 
energy,  E . Recent  developments  in  electronic  structure  theory  suggest 

that  many-body  perturbation  theory  provides  an  efficient  method  for 
computing  the  electronic  part  of  the  critical  energy.  Calculations 

for  the  formyl  radical  predict  a critical  energy  of  18.3  kcal  mol-1, 

which  agrees  well  with  the  value  assumed  in  this  research1'7.  A 

1 8 

forthcoming  report  will  discuss  these  calculations  in  detail 


if.  U.l;.  Adams,  G.H.  Bent,  G.D.  Purvis  and  R.J.  Bartlett,  manuscript 
in  preparation. 

18,  G.F.  Adams,  G.l>.  Bent,  G.D.  Purvis  and  R.J.  Bartlett,  ARKAIH'OM-TR 
in  preparation. 
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